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ABSTRACT 

We report the first detection of a sharp spectral feature in a Narrow-Line Seyfert 1 galaxy. Us- 
ing XMM-Newton we have observed IH 0707^95 and find a drop in flux by a factor of more 
than 2 at a rest-frame energy of ~7 keV without any detectable narrow Fe Ka line emission. 
The energy of this feature suggests a connection with the neutral iron K photoelectric edge, 
but the lack of any obvious absorption in the spectrum at lower energies makes the interpreta- 
tion challenging. We explore two alternative explanations for this unusual spectral feature: (i) 
partial covering absorption by clouds of neutral material and (ii) ionised disc reflection with 
lines and edges from different ionisation stages of iron blurred together by relativistic effects. 
We note that both models require an iron overabundance to explain the depth of the feature. 
The X-ray light curve shows strong and rapid variability, changing by a factor of four during 
the observation. The source displays modest spectral variability which is uncorrelated with 
flux. 
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1 INTRODUCTION 

IH 0707-495 (B = 15.7; z = 0.0411) has displayed some of the 
most rapid X-ray variability observed in a Seyfert galaxy (Turner 
, et al. 1999; Leighly 1999a). Its X-ray spectrum is steep and shows 
' a strong "soft excess" (Vaughan et al. 1999; Leighly 1999b) and its 
optical spectrum shows narrow H/3 (FWHM H/3 = 1050 km s"^) 
and strong Fe II emission (Remillard et al. 1986, Leighly 1999b), 
characteristic of many ultrasoft Seyferts. On the basis of FWHM 
H/3 IH 0707-495 is classified as a Narrow-Line Seyfert 1 (NLSl) 
galaxy. 

Such ultrasoft NLSls are of particular interest because their 
extreme properties are most likely produced by an extreme value 
of an underlying physical parameter (e.g. Pounds, Done & Osborne 
1995; Boiler, Brandt & Fink 1996; Brandt & Boiler 1998; Vaughan 
et al. 2001). However, with a few exceptions, their timing and spec- 
tral properties proved difficult to study in detail with previous X- 
ray observatories as a result of e.g., regular light curve interruptions 
from Earth occultation and limited spectral bandpass. 

IH 0707-495 was observed as part of a guaranteed time pro- 
gramme to study the timing and spectral properties of the X-ray 
brightest and most prominent ultrasoft NLS 1 s using XMM-Newton . 
In this Letter we report the first results from the EPIC data and 
in particular the detection of an unusual spectral feature at around 
7keV. 



2 OBSERVATION AND DATA REDUCTION 

XMM-Newton observed IH 0707^95 on 2000 October 21 during 
rev. 0159 for a duration of 46 ks, during which all instruments 
were operating nominally. The EPIC pn camera was operated in 
full-frame mode, and the two MOS cameras were in large-window 
mode; all three cameras used the medium filters. Extraction of sci- 
ence products from the Observation Data Files (ODFs) followed 
standard procedures using the XMM-Newton Science Analysis Sys- 
tem version 5.1.0 (SAS^. 

The raw MOS and pn data were processed to produce cal- 
ibrated event lists. Unwanted hot, dead or flickering pixels were 
removed, likewise events due to electronic noise, and event ener- 
gies were corrected for charge-transfer losses. The latest available 
calibration files^ were used in the processing. Light curves were 
extracted from these event lists to search for periods of background 
flaring and showed the background to be stable throughout the du- 
ration of the observation. The total amount of "good" exposure time 
selected was 42688 s for MOS 1, 42692 s for M0S2 and 38 106 s for 
the pn. (The pn exposure time includes a correction for out-of-time 
events in full-frame mode; Striider et al. 2001). 



* Seeptt p : / / xmm . vilspa. esa. es / u ser /sas_top .html 
^ http://xmm.vilspa.esa.es/ccf/ 
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Figure 1. Raw EPIC pn (blue) and MOS (red) count spectra of IH 0707^95 shifted to the rest-frame of the source. (The MOSl and MOS2 spectra were 
combined for display purposes only.) These spectra clearly show a significant drop at ~ 7 keV. The inset panel shows the pn spectrum, again in the rest-frame 
of the source, with the expected position of the neutral iron K edge mai'ked (dotted line). The bottom panel shows the data compared to a power-law model, 
with the slope fixed at F = 2.7. The model was not fitted to the data but chosen to emphasise the shape of the spectrum at both soft and hard energies. 



Source data were extracted from circular regions of radius 
60" for the MOS data, and radius 35" for the pn. (The smaller 
extraction region was used in the pn in order to avoid the CCD 
chip boundary ~ 40" from the pn aimpoint.) The total number of 
counts for the pn is 52377, for MOSl and M0S2 14501 and 14634, 
respectively. For the pn camera, both source and background re- 
gions were taken from inside the area corresponding to the hole 
in the pn circuit board. An examination of background images and 
spectra confirmed that the instrumental fluorescence radiation, es- 
pecially above 7 keV from Cu and Ni, is not affecting the source 
or background spectra. There are no signs of pile-up in the source 
region. Events corresponding to patterns 0-12 were extracted from 
the MOS data. Patterns 0-4 (single and double pixel events only) 
were used for the pn analysis. 



3 SPECTRAL ANALYSIS 

The source spectra were grouped such that each spectral bin con- 
tains at least 20 counts and were fitted using the XSPEC 11.0.1 
software package (Arnaud 1996). The latest publicly available re- 
sponses were used ^. The spectra from the three detectors were fit- 
ted simultaneously but with the relative normalisations free to vary. 
The quoted errors on the derived best-fitting model parameters cor- 
respond to a 90 per cent confidence level for one interesting param- 
eter (i.e. a Ax^ = 2.7 criterion) unless otherwise stated. Values of 

■t ml_medv9ql9t5r5_alL15.rsp and m2_medv9ql9t5r5_alL15.rsp for 
MOSl and MOS2. respectively and epn_ff20_sdY9_medium.rmf for the pn 



_ffo = 70 km s~^ Mpc~^ and go = 0.5 are assumed throughout, 
and fit parameters (specifically edge and line energies) are quoted 
for the rest frame of the source. 

Fig |l| shows the raw count spectra from the EPIC cameras. The 
broad-band spectrum is steep. A simple power-law fitted across the 
entire spectrum gives a slope V ~ 3.8, but this fit is dominated by 
the spectrum below 1 keV where the statistics are best. The most 
striking spectral feature is the sudden drop by a factor of more than 
2 at around 7 keV, visible in both pn and MOS spectra (note that this 
is a lower limit as the energy resolution of about 1 80 eV will smear 
a sharp edge). The width of the feature is not significantly larger 
than the energy resolution. The response matrices do not introduce 
any strong features around the observed energy of this feature, nor 
do the cosmic or instrumental backgrounds. In the 7-10 keV range 
the source is detected with a significance level of about 20 a, there- 
fore the spectral feature is intrinsic to the source. However, the lim- 
ited photon statistics above 10 keV do not allow us to constrain the 
expected recovery of the flux at higher energies. 

In the absence of any known systematic effect that could in- 
troduce such a feature into the spectrum, the rest of this paper dis- 
cusses possible physical origins of the feature assuming it to be 
intrinsic to the source spectrum. The much higher statistical weight 
at lower energies (due to the strong soft excess emission and peak 
in the EPIC response at ~ 1 keV) meant that fits to the full-band 
spectra were dominated by the soft emission. In order to concen- 
trate on the feature at 7 keV the data below 2 keV were ignored in 
the following analysis. 
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Figure 2. Spectral fits to the 2-10 keV EPIC data. The upper panels demon- 
strate the models (in Fi, units) and the lower panels show the data/model 
residuals. The left panels show the simple power-law and edge model, the 
middle panels show the partial covering model, and the right panels show 
the 3 X solar Fe ionised reflection model. 



3.1 Absorption by neutral iron? 

The most natural interpretation for a drop in the spectrum at just 
above 7 keV is as an absorption edge from a large column of neu- 
tral iron (7.112 keV, Henke et al. 1993; 7.111 keV, Bearden 1964). 
The 2-10 keV spectrum (see Fig |^ can be adequately parameter- 
ized {y^ — 192 /182 dof) using a power-law and a deep absorption 
edge (r = 1.8 ± 0.3). While providing an acceptable fit to the data, 
this model is unsatisfactory for a number of reasons. Firstly, the 
2-10 keV power-law slope is unusually flat (T = 1.07 ± 0.09). 
Secondly, the steep spectrum below 1 keV suggests there is little 
'cold' absorption above that expected from the Galactic column 
(A^H = 5.8 X 10^° cm"^ Dickey & Lockman 1990). Indeed, the 
column of cold iron expected from the depth of the putative edge is 
A^'pe « 6 X 10^^ cm"^ which corresponds to iVn « 2 x 10^" cm"^ 
for solar abundances. At this column density the spectrum should 
be dominated by absorption at low energies (L shell iron and K 
shell absorption from lower Z elements), whereas it actually steep- 
ens considerably in the soft band. Absorption by ionised material 
can produce a strong iron K edge without significant soft X-ray ab- 
sorption, but the energy of the putative edge suggests the absorber 
is not strongly ionised. The edge energy at 7.04 ± 0.07 keV is con- 
sistent with the neutral iron K edge energies given in the literature. 

The fluorescent yield for iron is close to 0.34 (Bambynek et 
al. 1972). However, there is no strong, narrow 6.4 keV emission 
line in the EPIC spectrum; the 90 per cent upper limit on the pho- 
ton flux in a narrow 6.4 keV line is 8 x 10~^ photons cm~^ s~^ 
(corresponding to an equivalent width of ~ 90 eV). The estimated 
photon flux absorbed by the putative edge is 1.2 x 10"^ photons 
cm~'^ s~^. For a spherically symmetric distribution of absorbing 
material one expects a ratio of fluorescence line to Fe-K absorbed 
flux approximately equal to the fluorescent yield, modified by the 
absorption of the 6.4 keV line (which is about a factor of 2.5 here). 
The measured ratio (of line/edge flux) is <_ 0.07, a factor ^ 2 
smaller than expected for a spherical distribution. This suggests 
that the column density along the line-of-sight is higher than the 
average surrounding the X-ray source. The solid angle subtended 
by the putative absorber, as seen from the X-ray source, is likely to 
be <j 2n sr (but see also section |5|). 

The following subsections discuss two alternative spectral 



models, specifically, partial absorption by neutral material and re- 
flection from ionised material. 



3.2 Partial covering model 

One possible explanation for the 2-10 keV spectral form of 
IH 0707^95 is partial covering (i.e., a patchy absorber; Holt et 
al. 1980): the observed spectrum is modified by cold absorption at 
high energies, while at lower energies a relatively unabsorbed com- 
ponent leaks through. A partial covering model does provide a good 
fit to the 2-10 keV spectrum (x^ = 177/178 dof), using three 
absorbed power-laws to explain the flat spectrum while the deep 
iron K edge in the most heavily absorbed component explains the 
7 keV feature (see Fig|^. A narrow 6.4 keV iron line is included in 
the model although, as noted above, it is only poorly constrained by 
the data. The best fitting model parameters are as follows: intrinsic 
power-law slope F — S.Slg ^, line-of-sight column densities (and 
covering fractions) of Nh ~ (0.04), 1.3 (0.12) and 4.6 (0.84) 
(in units of lO'^^ cm~^) and an iron abundance of about 35. Partial 
ionisation of metals could help, but the sharp spectral edge sug- 
gests that iron is not strongly ionised and iron L-shell absorption 
is unavoidable. The edge energy in the partial covering model is 
(7.10 ± 0.04) keV consistent with the energy of the neutral iron K 
edge. Fitting the 0.3-10 keV spectral energy distribution by adding 
a black body component, does not affect the spectral parameters 
derived above. 



3.3 Ionised reflection model 

A strong Fe edge may be found in the reflected emission from 
an irradiated accretion disc which is not too highly ionised (Ross, 
Fabian & Young 1999). To investigate if such models can describe 
these data, the constant-density reflection models of Ross & Fabian 
(1993) were fitted to the 2-10 keV spectrum. Since the reflected 
emission may originate close to the central black hole, relativistic 
blurring was applied to the model during fitting (using the LAOR 
code; Laor 1991). Reflection from ionised material was strongly 
preferred over neutral reflection (Ax^ ~ —40), an iron abundance 
approximately 3 times solar was needed (Ax^ = — 10 from a solar 
abundance model), and a Kerr space-time geometry was preferred 
over a Schwarzschild one (rin = 3^, rout= 8.3 gravitational radii). 
However, a model with a reflection fraction of unity can only give 
an adequate fit to the data (x^ ~ 258/180 dof) between 2 and 
10 keV, and it cannot fully account for the sharp spectral drop at 
~7 keV. Allowing the reflection fraction to increase to 9 resulted 
in a steeper continuum (F — 1.68 from 1.54) and an improved fit 
(X^ = 201/179 dof; see Fig||). The high value of R can be de- 
creased by increasing the iron abundance of the reflector. The ioni- 
sation parameter is ^ ~ 750 erg cm s~^. At this level of ionisation 
the reflection-dominated model predicts a strong He-like Fe Ka 
line with an equivalent width of 1.8 keV. This model clearly is not 
as good at describing the data as the partial covering model. It is 
extremely difficult for ionised reflection to account for the depth of 
the drop at ~7 keV without invoking a very extreme Fe abundance 
and/or reflection fraction. 



3.4 Other models 

The blue wing of a very strong relativistically broadened and highly 
redshifted iron Ka line can produce a drop at 7 keV. A power- 
law plus LAOR line (rest energy 6.4 keV) provides a good fit 
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Figure 3. EPIC pii light curve in the 0.1-10 keV energy range with a bin 
size of 200 s, demonstrating the strong and rapid variability. The "high flux" 
and "low flux" intervals (see section U) are marked differently. 



Figure 5. Hardness ratios using 1000 s time bins. The first two hardness 
ratios (at the top) are acceptable fits to a constant. The third is unacceptable 
at a ~ 95 per cent level, the last one is unacceptable at > 99.9 per cent. 
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Figure 4. RMS variability versus energy for the EPIC pn light curve. 
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Figure 6. Ratio of "low flux" and "high flux" spectra. 



(X^ = 182/180 dof) with an inner disc radius r-m < 2rg and an 
underlying power-law slope of F = 2.0815)^13. However, the very 
high equivalent width (EW — 5 keV) required to explain the size 
of the 7 keV drop is difficult to justify. Alternatively, a reflection 
spectrum from cold material will contain a neutral iron edge. Fit- 
ting with a PEXRAV model gives a good fit (x^ = 192/182 dof) 
with a reflection fraction > 200 (essentially just reflection) and an 
iron overabundance > 3. The dominance of the reflected compared 
to the primary emission, and the weakness of the associated iron 
Ka line appear fatal for this interpretation. 



4 TEMPORAL ANALYSIS 

Fig ^ shows the 0.1-10 keV EPIC pn light curve in 200 s bins 
(The MOS light curves are essentially identical to the pn light 
curve). The source again showed strong (factor of ~4 change dur- 
ing the observation) and rapid variability, as previously seen by 
ASCA (Leighly 1999a). The relatively low apparent luminosity 
of the source during this observation, I/0.2-10 ~ 8 x 10^^ erg 
s~^, means that the rapid variability translates to a fairly modest 
rate of change of luminosity. For example, the rapid rise toward 
the end of the observation corresponds to a luminosity increase of 
AL/ At ~ 2.4 X 10"^^ erg s~^. The corresponding radiative effi- 
ciency, using the argument of Fabian (1979), is only rj > 0.2 per 
cent. 

Light curves in various energy bands were analysed to search 
for spectral variability. The variations in each band appear very 
similar. The fractional variability amplitude Fvar (see Edelson et 
al. 2001) was measured in different energy bands (Fig ^ using 
light curves binned to 1000 s and is consistent with a constant 
(fractional) variability amplitude. A cross-correlation analysis us- 



ing the discrete correlation function (DCF; Edelson & Krolik 1988) 
showed the light curves to be correlated at zero-lag, with no evi- 
dence for any time lags. 

Hardness ratios were also examined to search for spectral vari- 
ability (see Fig These do show spectral variability (the lower 
two panels) but are uncorrelated with flux. As a final check of 
(flux-correlated) spectral variability separate pn spectra were ex- 
tracted from time intervals when the source flux was below the 
mean ("low flux" spectrum) and above the mean ("high flux" spec- 
trum). The ratio of the two spectra (Fig ^ was consistent with a 
constant, re-enforcing the claim of no flux-dependent spectral vari- 
ability. Unfortunately, the limited number of counts above 7 keV 
(~ 140 source counts in the pn) mitigates against a detailed analy- 
sis of the variability of the 7 keV spectral feature. 



5 DISCUSSION 

The XMM-Newton spectrum of IH 0707^95 has revealed a sharp 
drop, by a factor ^ 2 in the spectrum at an energy just above 7 keV. 
This is to our knowledge the first detection of such a feature in an 
AGN. The limited statistics and complex continuum mean that the 
identification of this feature remains ambiguous. We have consid- 
ered models involving either partial covering of the X-ray source 
by cold material, or emission from an ionised accretion disc, but 
both these explanations have drawbacks. 

5.1 Comparison of the two models 

The energy of the feature is strongly suggestive of absorption by 
neutral iron. However, a complex partial covering model is required 
in order that the soft X-ray spectrum be relatively unaffected by the 



© 2001 RAS, MNRAS 000, 



XMM-Newton discovery of a sharp spectral feature in IH 0707-495 5 



absorber. A further implication of tlie partial covering model is that 
the underlying power-law is unusually steep (F ~ 3.5). A similar 
feature, namely a sharp drop near 7.1 keV with relatively little flu- 
orescence emission, was observed in the peculiar X-ray binary Cir 
X-1 (Brandt et al. 1996). In this case the spectral feature is almost 
certainly due to partial covering. 

An important consequence of the partial covering model is 
that the unabsorbed X-ray luminosity of IH 0707-495 is an order 
of magnitude higher than the observed (absorbed) luminosity. We 
have calculated the Oox values for the absorbed and unabsorbed 
fluxes at 2 keV (the flux density at 2500 A is taken from Leighly 
2000). The values are Qox = —2.0 (absorbed X-ray flux) and 
Oox = —1.6 (unabsorbed). This is suggestive of absorption (see 
Brandt, Laor & Wills 2000; Brandt et al. 2001, Figure 3) as almost 
all the objects with Qox < —1.8 suffer from heavy UV and X-ray 
absorption. However, there is no evidence for intrinsic absorption 
in the HST spectrum of IH 0707-495 (Leighly 2000). 

The ionised disc model including relativistic blurring seems to 
provide another possibility for the sharp spectral feature. The blue 
horn of the relativistically broadened Ka line emission, combined 
with the iron edge can produce a drop around 7 keV. However, there 
still remain some residuals above 6.8 keV in the observers frame 
and it appears somewhat surprising that the energy of the feature is 
almost exactly that of the neutral iron edge. 



5.2 A physical scenario for the partial covering model 

The simplest configuration for the absorber would involve the ab- 
sorber completely covering our line of sight, the unabsorbed flux 
being a scattered component from some other direction. This scat- 
tered component would not then vary rapidly and so would cause 
spectral variability when combined with the (variable) absorbed 
flux. If on the other hand the absorber does partially cover the 
source then it presumably is at most only a few 100 light-seconds in 
size and its density exceeds ~ 10^^ cm~^. Photoionization consid- 
erations require that the cloud thickness be less than ~ lO^^rij cm, 
where the distance of a cloud from the nucleus is rislO'^^ cm. The 
source position, size and spectrum must remain fixed while the flux 
is rapidly varying in order to avoid strong spectral variability. 

These problems are alleviated if the absorbers are close to, or 
even within, the emission region. The absorbing clouds must then 
be very dense to remain cool and therefore held by magnetic fields. 
The possible presence of such clouds in AGN in general has been 
discussed by Rees (1987); Celotti, Fabian & Rees (1992); Kun- 
cic, Blackman & Rees (1996); Kuncic, Celotti & Rees (1997) and 
Malzac (2001). Brandt & Gallagher (2000) first suggested that such 
clouds might be particularly relevant for NLSls. Cool gas trapped 
by the magnetic field is compressed to extreme densities by the 
high radiation pressure. Clouds and filaments of such cold gas may 
accompany the active regions above an accretion disc. 

This possibility has been investigated using a relativistically- 
blurred model involving emission and partial covering. In order that 
the 7.1 keV feature remains sharp in the model, the inclination must 
be <. 15° . If the absorbing material lies outside the emission region 
(say in a toroidal structure at 10^^ cm) then only the line-of-sight 
material affects the spectrum and this constraint is relaxed. 

It is interesting to speculate whether the absorbing material is 
connected with the radiation-driven outflow suggested by Leighly 
(2000) in IH 0707-495, and whether it affects the observed rapid 
X-ray variability. Clearly it is crucial to search for similar fea- 
tures in other objects. The XMM-Newton spectrum of PG 121 1-1-143 



(Reeves et al. in prep.) appears to show a similar drop at ~ 7.3 keV. 
If such features are common in other objects, this will further con- 
strain the range of suitable models. 
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